Entangling the free motion of a particle pair: an experimental scenario 
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The concept of dissociation-time entanglement provides a means of manifesting non-classical cor- 
relations in the motional state of two counter-propagating atoms. In this article, we discuss in detail 
the requirements for a specific experimental implementation, which is based on the Feshbach disso- 
ciation of a molecular Bose-Einstein condensate of fermionic lithium. A sequence of two magnetic 
field pulses serves to delocalize both of the dissociation products into a superposition of consecutive 
wave packets, which are separated by a macroscopic distance. This allows to address them sep- 
arately in a switched Mach-Zehnder configuration, permitting to conduct a Bell experiment with 
simple position measurements. We analyze the expected form of the two-particle wave function 
in a concrete experimental setup that uses lasers as atom guides. Assuming viable experimental 
parameters the setup is shown to be capable of violating a Bell inequality. 
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I. INTRODUCTION 

The increasing control of quantum systems that has 
been obtained in the last decades brings into reach the 
possibility of performing hitherto unprecedented tests of 
the principles of quantum mechanics In partic- 

ular, there is a growing interest in demonstrating non- 
classical behavior in systems or with respect to features 
that would be assigned to the classical realm accord- 
ing to our every-day experience. One exciting approach 
aims at luring mesoscopic mechanical oscillators into the 
quantum regime, see [6|] and refs therein. Another ap- 
proach focuses on the generation of highly non-classical 
states of continuous variables with direct classical cor- 
respondence, such as creating superpositions of macro- 
scopically different phase space coordinates ^-'T^. The 
recent achievements in the generation and manipulation 
of ultracold atoms and molecules (ll| - [l6j suggest to im- 
print non-classical features even on the free motion of 
material particles. For example, it has been proposed 
to perform the original Einstein- Podolsky- Rosen thought 
experiment based on the dissociation of a molecular Bose- 
Einstein condensate (BEC) [13, El- 

In a recent article we go one step further by proposing 
to violate a Bell inequality in the motional state of ma- 
terial particles using position measurements in the end 
[l9| . The two- particle state considered in this proposal 
can be characterized by its generation procedure. Start- 
ing from a diatomic Feshbach molecule, which is trapped 
in an optical dipole trap, a sequence of two short mag- 
netic field pulses creates the superposition of two subse- 
quent dissociations. Due to momentum conservation the 
dissociated atoms propagate in opposite directions along 
some wave guide, with each atom delocalized into a pair 
of macroscopically separated, consecutive wave packets, 
corresponding to the early and the late dissociation in- 
stant, see Fig. [TJ An observer cannot tell whether the 
dissociation took place at the early or at the late instant, 
while both atoms must have started their journey at the 



same time. In this sense, the generated two-particle state 
may be called dissociation-time entangled (DTE). 

Since each atom in a DTE state is described by two 
spatially distinct propagating wave packets, which corre- 
spond to the early and the late dissociation times, the use 
of switched, asymmetric Mach-Zehnder interferometers 
can serve first to endow these wave packets with addi- 
tional relative phases and then to rejoin them. In detail, 
the early wave packets propagate through the long arms 
of the interferometers, whereas the late ones are directed 
through the short arms, see Fig. [T]c). If the detour of 
the early wave packets is chosen appropriately, they re- 
join the late ones in the final beam splitters. Provided the 
two dissociation processes are coherent in the sense that 
no information about the dissociation time has leaked to 
the environment, the rejoined wave packets may inter- 
fere. The joint probability for detecting the particles in 
the output ports of the two interferometers then exhibits 
a fringe pattern under variation of either path length dif- 
ference, which is in principle capable of violating a Bell 
inequality. We emphasize that the only relevant informa- 
tion required to extract these correlations is the output 
port combination where the atoms are detected. This 
information can be obtained by simple position measure- 
ments, without the need of a particular spatial or tempo- 
ral resolution. And since the position measurements can 
be understood from a classical physicist's point of view, 
they are indeed in line with the above mentioned pursuit 
of focusing on features that have a classical analogoue. 

We note that a similar Bell test procedure was first 
introduced for photons in [l^, HH . This led to the ex- 
perimental demonstration of a Bell inequality violation 
for parties separated by several kilometers, albeit using 
a postselection procedure required due to the lack of suf- 
ficiently good switches for photons. The implications for 
material particles, in particular the detrimental effect of 
dispersion, are discussed in detail in p^ . 

In this article we devise a concrete scenario to perform 
the described DTE Bell test, which is sufficiently detailed 
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FIG. 1: (color online) a) Experimental setup for the DTE Bell test. It starts with a ^Li2 Bose-Einstein condensate of about 10^ 
molecules in an optical dipole trap. A sequence of two short magnetic field pulses, as shown in b), dissociates on average one 
molecule per shot into a pair of counter-propagating atoms. The resulting two-particle state corresponds to a superposition of 
two subsequent dissociations, where each atom is delocalized into two consecutive wave packets that are separated by a distance 
of 5mm. The two atoms propagate in opposite directions along the laser guide and eventually enter switched, asymmetric Mach- 
Zehnder interferometers that are implemented by two more laser beams crossing the guide in a triangular arrangement at small 
angles. These second lasers are switched on after the early wave packets have passed, but before the late wave packets arrive. 
An additional perpendicular, blue-detuned blocking laser beam, which is not depicted in the figure, may be used to support 
the deflection of the late wave packets. This way, the early wave packets have to take the long paths, whereas the late ones 
are directed into the short paths, as depicted in c). If the path length differences cancel the distance between the early and 
the late wave packets, they interfere and the joint probability of detection in the output ports of the interferometers exhibits 
a fringe pattern as a function of the path length difference variation. Correlations of the single-particle detection events thus 
have the potential to violate a Bell inequality. 



and realistic to suggest that violating a Bell inequality 
should be possible within reach of present day technol- 
ogy. Section |TT] provides a detailed quantitative account 
of the proposed setup, based on realistic choices for the 
experimental parameters. It involves pulse separations 
on the order of seconds and implies atom velocities on the 
order of centimeters per second, as well as a truly macro- 
scopic spatial separation of the single atom wave packets 
on the order of centimeters. In Section IIIII we discuss 
and numerically evaluate the two-particle superposition 
state that is generated in the wave guide in the course of 
a double dissociation with two square pulses. Section ITVl 
demonstrates that this state, if exposed to the switched, 
asymmetric Mach-Zehnder processing, generates a fringe 
pattern which violates a Bell inequality. Finally, we con- 
clude in Section |V] by pointing out the advantages of the 
proposed scheme. 

II. EXPERIMENTAL SCENARIO FOR A DTE 
BELL TEST 

In the following we discuss quantitative details of a 
possible experimental setup that should be capable of 
generating a pair of dissociation-time entangled atoms 
and to violate a Bell inequality in terms of single-particle 



interferometry and simple position measurements. We 
work at this example to prove the viability of the pro- 
posal presented in [T9| . but of course many alternative 
implementations are conceivable. We emphasize that the 
setup is designed to allow for time and length scales on 
the order of seconds and centimeters, establishing entan- 
glement with respect to single-particle properties that are 
truly macroscopically distinct. In this section we focus 
on a consistent set of experimental parameters, deferring 
the theoretical analysis of the parts concerning the dis- 
sociation process to Section IIIII 

We suggest to use a dilute molecular Bose-Einstein 
condensate (BEC) produced from a 50:50 spin mixture of 
fermionic ^ Li as a starting point. Such a fermionic mix- 
ture seems favorable compared to bosonic ingredients, 
since huge lifetimes of more than 10s can be achieved 
due to Pauli blocking of detrimental 3-body collisions 
p^ . This makes truly macroscopic time separations t 
between the two dissociation pulses conceivable, and we 
choose T = Is in the following. It has been demon- 
strated that the molecular ^Li BEC can be prepared effi- 
ciently and with near-perfect purity [l2] , and the compar- 
atively small mass of lithium reconciles reasonable prop- 
agation velocities, on the order of Icm/s, with resolvable 
de Broglie wave lengths, on the order of 10/im. 

The BEC is prepared in an optical dipole trap created 
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by two red-detuned, far ofF-resonant, crossing laser beams 
with a wave length of, say A = l^m, see Fig. [T]a). The 
guiding laser, which connects the two single-atom inter- 
ferometers, serves as an atom guide which compensates 
the gravitational force. Moreover, it must be chosen suf- 
ficiently strong that the transverse motion of the atoms 
remains frozen in the ground state. This way an effec- 
tively one-dimensional description of the atomic motion 
can be used. 

The following numerical values are chosen to be con- 
sistent with a pulse separation of r = Is and a rela- 
tive velocity of Wroi = Icm/s. We assume that the guid- 
ing laser is given by a Gaussian laser beam and oper- 
ated at a power of Pq — 32.85W with a (center) waist 
of lUo.G = 216/im, such that the optical dipole poten- 
tial resulting from a molecular polarizability of aLia (A = 
1/im) = 88. SA has a (center) trap depth of 30/iK. It 
thus is strong enough to compensate gravitation, sup- 
porting several thousand transversally bound states. The 
implied Rayleigh length of about Zq.g ~ 15cm sets the 
scale for the extension of the whole setup. Decoherence 
due to photon scattering can be estimated to occur at a 
rate on the order of 0.05s^^, which limits the overall du- 
ration of a single experiment to remain below 20 s. The 
implied transverse trap frequency of ujq /2Tr — 300Hz, on 
the other hand, limits the maximum dissociation veloc- 
ity not to exceed about Urei = Icm/s, if we demand the 
transverse motion of the atoms to remain frozen in the 
ground state. The corresponding de Broglie wave length 
Arei = 12.4/im, which sets the scale for the fringe separa- 
tion, is in compliance with viable stability requirements 
for the interferometers. 

The second laser beam intersects the guiding laser per- 
pendicularly, creating an elongated dipole trap for the 
BEG within the laser guide. At the same time, it must be 
sufficiently shallow to be overcome with reasonable mag- 
netic field pulses. Moreover, the resulting longitudinal 
trap ground state IV't), which is supposed to determine 
the center of mass state of the trapped molecules, must 
meet stringent conditions regarding its momentum spec- 
trum. Only a sufficiently narrow momentum spectrum 
guarantees that the dispersion-induced reduction of the 
non-local correlations due to the momentum spread in 
the center of mass motion admits the violation of a Bell 
inequality [l^ . This demands a comparatively small trap 
frequency. Specifically, for an atomic propagation veloc- 
ity of 5mm/s, the trap frequency should not greatly ex- 
ceed WT/27r = 0.25Hz. 

A reasonable trap depth of UT/k-Q = 50nK combined 
with the required trap frequency of ujt/2ti — 0.25Hz im- 
plies the trap laser (center) beam waist to be ?«o,t = 
1.1cm. By employing an elliptic Gaussian laser beam 
with waist ratio Wx^r/wy^T — 10, the required laser power 
can be kept reasonably low, at Pt = 13. IW. The corre- 
sponding photon scattering rate of 10~'*s~^, on the other 
hand, does not impose any limitations. 

The preparation is arranged such that only a small 
number of molecules, on the order of 10^, remains in 



the BEG at the end. These can be taken to be non- 
interacting, so that the initial longitudinal center of mass 
state of the molecules indeed is given by the trap laser 
ground state, parametrized by the trap frequency wt- 

A DTE atom pair is extracted from the BEG by apply- 
ing a sequence of two short square pulses in the external 
magnetic field, as shown in Fig. [T]b). The pulses must 
be chosen such that the desired dissociation velocity of 
Vj:c\ ~ Icm/s is accomplished. To this end the pulses not 
only have to provide the corresponding kinetic energy 
but also the energy to overcome the longitudinal trap 
depth of Ut = fceSOnK. Moreover, the generated state 
of the relative motion must meet similarly stringent con- 
ditions concerning its momentum spectrum as the trap 
laser ground state IV't)- 

The desired dissociation velocity of v^ei = Icm/s is ob- 
tained for a narrow, sufficiently isolated resonance with, 
e.g., resonance width Ai?rcs = ImG and /iics = O.Ol/ie, 
by applying a pulse duration of T = 60ms and a pulse 
height of Bq + AB — B^cs = 200mG. Here, ^res denotes 
the difference between the magnetic moments of the reso- 
nance state and the open channel, Bq stands for the mag- 
netic field base value before and after the pulses, and i?rcs 
indicates the resonance position. In particular, the gen- 
erated state meets the tight restrictions on the momen- 
tum distribution as imposed by dispersion. In Section Hill 
we will analyze the momentum spectrum resulting from 
such a double pulse sequence and explicitly demonstrate 
its capability of violating a Bell inequality. 

The dissociation probability per double pulse sequence 
depends, apart from the shape of the magnetic field 
pulses, on the resonance parameters ABj-es and /ires, on 
the guiding laser trapping frequency ujq , and on the back- 
ground channel scattering length abg. For a generic value 
of Obg — lOOoo the dissociation probability amounts to 
a few percent. We thus rely on post-selection in order 
to guarantee that only one atom pair per shot enters the 
interferometers. This is however unproblematic, since 
the final fiuorescence detection of the slow, strongly con- 
fined atoms can be done with single particle resolution, 
so that it is easy to disregard the cases of too many atom 
pairs in the process. For a start, we stick to this sim- 
ple post-selection procedure; in a more refined setup it 
is conceivable to use a specially prepared optical lattice 
where each site is occupied by at most one molecule [l^ . 

After the completion of the dissociation sequence, and 
for our chosen time separation of r = Is between the 
early and the late dissociation pulse, the corresponding 
early and late wave packets of each particle propagate at 
a velocity of 5 mm/s, separated by a distance of 5 mm on 
each side. This constitutes a truly macroscopic delocal- 
ization of each atom. Immediately after the dissociation 
process, the widths of the early and late single-particle 
wave packets are on the order of about 200/im. Their 
narrow momentum spectra guarantee that these wave 
packet extensions are not appreciably modified during 
the propagation to the interferometers if the propaga- 
tion time does not exceed about 10s. The early and late 
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wave packets are thus spatially still sufSciently distinct 
when arriving at the interferometers for the switches to 
be applicable. 

The interferometers are implemented by two more red- 
detuned laser beams crossing the guide in a triangular ar- 
rangement at small angles, see Fig. [TJ The path length 
differences must compensate the 5mm distance between 
the early and late wave packets of each particle. Small 
variations of the path length differences on the order of 
the de Broglie wave length Aroi = 12.4/im result in a vari- 
ation of the relative phase between the early and the late 
wave packets, as required for the Bell test [l^. The ad- 
ditional laser beams are switched on only after the early 
wave packets have passed the crossing points, but before 
the late ones arrive there. The crossings of the beams 
then act as beam splitters, while the required atom mir- 
rors may be realized using evanescent light fields or blue- 
detuned laser beams perpendicular to the interferometer 
planes (24,, 25.]- Such perpendicular blocking beams at 
the crossing point of the first beam splitters could also 
support the full deflection of the late wave packets into 
the short arms. This way, the early wave packets pass 
through the long arms, whereas the late wave packets 
are deflected into the short arms, as depicted in Fig. [1] 
c). Note that a simplified setup could do without the 
switching, replacing the switches by ordinary beam split- 
ters, at the cost of 50% post-selection. 

Resonant laser beams crossing the guiding lasers be- 
hind the final beam splitters, as sketched in Fig. [T]a), 
effect fluorescence detection with near unit efficiency and 
with single particle resolution. No particular spatial or 
temporal resolution of the detection is required, since the 
only relevant information for the Bell test is the combi- 
nation of beam splitter output ports where the particles 
are detected. Similar to the original Bell test based on 
spins, the joint probability distribution for detecting the 
particles in particular output port combinations shows an 
interference pattern under variation of either path length 
difference, which can violate a Bell inequality. Decoher- 
ence due to photon scattering in the guiding laser beam 
sets the scale for the distance between the BEC trap and 
the detection lasers to be on the order of 10 cm in our 
example. 



III. DTE STATE FROM A DOUBLE 
DISSOCIATION PULSE 



persion. 

A two-channel single-resonance calculation shows that 
after an arbitrary magnetic field pulse sequence (close to 
an isolated resonance) the dissociated part of the state, 
\^hg{t)), is described, for low energies, at positions far 
from the dissociation center, and at large times, by the 
asymptotic form 



(0), 



Here j denotes the free (two-particle) propagator in 
the longitudinal direction, while the transverse motion 
is frozen in the harmonic ground state, It^o™) ^^-d |</3"o)' 
resp., of the guiding laser beam. The longitudinal part of 
the two-particle state, expressed in the basis of the center 
of mass and relative momenta Pcm and Pid, is determined 
by [19,] 



(Pcin,Prcl|W^) = -5- (PcmlV'T), 

(1) 

where C{uj) is the Fourier transform of the closed channel 
probability amplitude C(i), (5(0;) = /_!°^ die'"*C(i). The 
quantity 
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normalizes the spectrum. The dissociation probability 
can be estimated by 
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(3) 



which involves the background scattering length Obg, the 
resonance width AB^os, and /ires, the difference between 
the magnetic moments of the resonance state and the 
open channel. As mentioned in Section we presume a 
narrow resonance width of AB^^s = ImG and short pulse 
durations T — 60ms, such that only a small fraction of 
the condensate dissociates. The coupling between the 
two channels is then sufficiently weak so that the back- 
action from the open channel to the bare resonance state 
can be safely neglected. In this case, a time-dependent 
magnetic field B{t) determines the closed channel prob- 
ability amplitude according to 



We proceed to derive the asymptotic form of the two- 
particle state that is generated by a double pulse se- 
quence in the described setup, and we verify that it has 
the DTE structure as described in Section H] In particu- 
lar, we will be able to specify the corresponding momen- 
tum spectrum, which plays a crucial role for the feasi- 
bility of the Bell test. Below in Section IIVI we will use 
this result to demonstrate explicitly that the experimen- 
tal scenario presented above yields a DTE state which 
can violate a Bell inequality in spite of wave packet dis- 



C{t) = C{to) exp 



dt'[E,,s{Bit'))-2UT + hujG] 



(4) 

The resonance state energy i^ics (-6(0 ) = Mics(S(t)— Bros) 
is measured with respect to the background channel con- 
tinuum threshold, whereas the trap laser potential depth 
Ut and the transverse trap frequency of the guiding laser 
loq cause an additional energy offset that has to be over- 
come in the dissociation process. C(to) indicates the 
probability amplitude for the two particles to be found in 
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the resonance state at some initial time to well before the resonance Bros, the corresponding background channel 

dissociation pulses take place. Since we can presume that component is negligible and we can safely set C(to) — 1 

the magnetic field before to has remained well below the from now on. 

I 



For a sequence of two square pulses, each of height AB and duration T, separated by the time r, 
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the Fourier transform of @ can be evaluated analytically, yielding 

T/irosABsinc [(w - /Xres[Bo + AB- B,,,]/h + 2UT/h - loq) T/2] 



?IL0 - Hrcs[Bo - -Bros] + 2Ut - tu^G 
j[2;7TT-/i„BABT+/^jca(-B„B-So)T-?itJG 



(5) 



(6) 



r]/n| 



r 



where we have neglected an irrelevant global phase. It is 
now helpful to introduce the mean energy 



/iros(Bo + AB- Bros) - 2L/t - hujQ 



the characteristic width 

= 2mh/T, 

and the pulse energy 

p^/m = fircsAB. 
The dissociated state H]) thus takes the form 

||C||(pcm,Prel|*.) = Tp^ipcJi^T) (7) 

., sinc[fe/4 + pg^i~pg) /Ap^] 

where the relative phase between the early and the late 
dissociation component is given by 

(j)r = [2UtT - flrcsABT + ^res(Brcs - Bo)T]/h + UqT. 

For the choice of parameters presented in Section [IT] 
the momentum spectrum ([7]) is sharply peaked at 
(PcmjProi) = (0, ±po), with po / m = 5mm/s. This 
is demonstrated by considering the ratios (Tp x/Po a-nd 
Ap/po, where (Tp,T denotes the momentum uncertainty 
exhibited by the longitudinal trap ground state |V't), 
reading in the harmonic approximation dp^T = y/huTin. 
This yields the ratio ap^^/po = 0.024, while the pulse 
duration T — 60ms implies Ap/po = 0.11. We can ex- 
ploit the smallness of these ratios in order to calculate 
the normalization ||C|p according to For this, we 
can focus on a single-pulse contribution, since in the 
limit T ^ T the spectrum ([T]) represents the superpo- 
sition of two spatially distinct dissociation states, which 



contribute equally to the norm. Within the range Up^T 
where (pcmlV'T) is non-negligible, the Feshbach contri- 
bution C{pl^/4:mh + p'^^i/mh+ 2ujq) to the right hand 
side of ([7]) is only weakly dependent on the center-of- 
mass momentum Pcm, which permits to replace pcm by 
its mean value pcm = 0. The remaining integral over 
the relative momentum proi can be evaluated in good ap- 
proximation by linearizing its argument in the vicinity of 
Po- Setting Prci ~Po ~ 2po(Prci -Po) within the region of 
non-negligible C{p'^^^/mh+ 2ujq) thus yields 



|C||2 = 27rr2ApVP0. 



(8) 



With this, and realizing that the momentum dependent 
phase factor on the right hand side of ([7]) effects a free 
time evolution, we end up with a longitudinal asymptotic 
wave packet of the form 



U 



(0) 



(9) 



Clearly, it is the superposition of an early and a late 
dissociation contribution. Note that the additional free 
time evolution of the early state component effects a 
dispersion-induced distortion with respect to the late 
state component. The inferred single-pulse momentum 
spectrum reads 



(Pcm,Prcl|*o) = 



Pop smc 



pD 



^ApiplJA^ 

x(pcm|^T>- 



Pi 



(10) 



The state \^o) describes a pair of symmetrically counter- 
propagating particles, in the sense that each particle can 
propagate in both directions, provided the other particle 
takes the opposite direction. The corresponding momen- 
tum distribution |(pcm,Prei|*o)P is plotted in Fig. [5] for 
our choice of parameters. As anticipated, the momentum 
spectrum ((TU)) is sharply peaked at (pcm,Prei) = (0, ±Po)- 
This is of great importance for the viability of the DTE 
Bell test with material particles, as was shown in ^19.] . 
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FIG. 2: (color online) Momentum distribution j (Pcm, Preil^'o) P for a single Feshbach dissociation pulse with pulse duration 
T = 60ms, magnetic field base value Bros — Bo ~ 200mG and pulse height AB — 400mG. Assuming a longitudinal trap depth of 
Ut /kB = 50nK and a transverse trapping frequency of u)g/2h = 300Hz this yields the dissociation velocity po/m = 5mm/s. The 
upper plot shows the momentum distribution in the pcm = plane on a logarithmic scale. As required, the distribution exhibits 
sharp peaks at (pcm,Prci) = (0, ±Po), corresponding to symmetrically counterpropagating particles. The lower logarithmic 
contour plot focuses on the vicinity of the (pcm, Prci) = (0, +po) peak, including the dependence on the center of mass momentum 
Pcm. The width in p^oi is characterized by Ap^ = 2mh/T, with (Ap/po)^ = 0.012, whereas the width in pcm is essentially 
determined by the momentum uncertainty cTp.T of the longitudinal trap laser ground state \iI)t}, which follows from the trap 
frequency ujt/2-k = 0.25Hz, yielding ap^T/po ~ 0.024. 



The double dissociation pulse sequence generates a 
two-particle state of the form where due to the ad- 
ditional free time evolution of the early state component 
the corresponding wave packets trail each other. The 
early and late wave packets are sufficiently spatially sep- 
arated such that they can be subjected to the required 
interferometric transformation. The double pulse disso- 
ciation state (O thus has the structure of a DTE state as 
described in Section HI The dissociation probability for 
the double pulse sequence follows from Eqs. ^ and 
For our parameters, we get |CbgP ~ 0.04, which confirms 
the premised small dissociation fraction and thus justi- 
fies a posteriori the weak coupling approximation in the 
derivation of (0]). 



mation described in Section U produces an experimen- 
tally resolvable fringe pattern that violates a Bell inequal- 
ity. In brief, the basic idea behind the interferometers is 
to endow the early state components with an additional 
phase before recombining them with the late state com- 
ponents. The feasibility of this procedure requires the 
early and the late wave packets not to overlap before ar- 
riving at the first, switchable mirror. If this condition is 
met by the DTE state ^ the joint probability for de- 
tecting a pair of counterpropagating atoms in a particu- 
lar output port combination labeled by (Ti, (72 {o'i = ±1) 
follows from 



IV. DTE BELL TEST 



In this section we demonstrate explicitly that the DTE 
state (O, when subjected to the interferometric transfor- 
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1 + ai(j2Re < e ' 



dprel|(Pcm,Prcl|*o'^^)P 



(11) 



X exp /i^£i^i(fl±M + iMlfl^ _ . T(pg^/4 + pg^i) 
2?i /i mh 



as shown in [19|. The two-particle momentum distribu- 
tion I (pcm J Prel 

der ives from the symmetric dis- 
tribution I (pcmiPrcil^o) P by restricting each particle to 
propagate into a particular direction, say, particle 1 into 
positive direction and particle 2 into negative direction. 
Hence we have 

|(Pcm,Prcl|*o"^^)P = 26l(p„l)|(pcm,Prcl|1'o)P, 

where the additional factor 2 is imposed by normal- 
ization. If the two-particle momentum distribution 
|(Pcm,Pici|^o)P is sufficiently well-behaved and the path 
length differences £i, £2 are chosen such that the early 
and late wave packets overlap in the output ports, the 
detection probability (jlip reveals a fringe pattern under 
variation of £1 and £2 ■ It is remarkable that the quality of 
this fringe pattern is only affected by the free time evolu- 
tion between the early and the late dissociation, whereas 
any subsequent time evolution does not have any effect, 
despite of the ongoing broadening of the wave packets. 

For the single pulse momentum spectrum (jlOp of the 
DTE state Q the joint detection probability pT|) can 
be evaluated numerically. The resulting fringe pattern 
for the output port choice cti = -1-1 and (T2 = +1 under 
variation of £1 can be seen in Fig. [3l 

Specifically, we choose £1 = £0 + A£ and £2 = —£o: 
where £0 = rpo/m defines the optimum overlap path 
length difference cancelling the separation between the 
early and the late wave packets. The variation of £2 or 
combinations of £1 and £2 would of course yield similar 
fringe patterns. The (constant) relative phase (j)r is ne- 
glected since it only results in a phase shift of the fringe 
pattern, while the upper and lower envelope remain un- 
affected. 

Inspecting Fig. [3l we can identify several generic fea- 
tures of the DTE fringe pattern. First, the characteristic 
length for its period is set by the de Broglie wave length 
of the relative motion Arei = h/po. This is natural, since 
we implement the relative phase shift by displacing the 
early and late wave packets with respect to each other, 
while the phase variation of a propagating wave packet 
is set by its de Broglie wave length. In our case, we get 
Aroi — 12.4/im. This means that, in order to be able to 
detect the fringe pattern, the interferometers have to be 
kept stable by about 1/ini. 

The envelope size of the fringe pattern, on the other 
hand, can be traced back to the unavoidable offset from 
the optimum wave packet overlap when varying the path 
length differences. One expects no interference if there 




FIG. 3: (color online) Joint probability Pdto(cri, cr2|^i, ^2) for 
detecting a pair of counter-propagating atoms in the output 
port combination cti = -1-1 and a2 — +1 under variation of the 
path length difference £1. At the outset, £1 and £2 are chosen 
such that optimum overlap of the early and the late wave 
packets is achieved in the output ports of the beam splitters. 
The upper plot shows the resulting fringe pattern as it would 
be seen by the experimenter. The characteristic length for 
its period is set by the de Broglie wave length of the relative 
motion, Aroi ~ 12.4/im. The finite range of the envelopes 
reflects the decreasing wave packet overlap with increasing 
offset A£ from the optimum overlap value. The lower plot 
highlights the corresponding fringe visibility. The dashed line 
marks the visibility threshold that has to be exceeded in order 
to violate a Bell inequality. A sufficient number of fringes 
exceeds the threshold to manifest the corresponding violation. 



is no overlap in the beam splitter output ports, result- 
ing in the uncorrelated probability distribution of 25% 
per output port combination. Indeed, the width of the 
fringe pattern, which is of about 200/im, matches approx- 
imately the initial position spread of the single-particle 
wave packets. Minor modifications to this can be traced 
back to the dispersion-induced distortion between the 
early and the late wave packets. 

However, the overall visibility reduction, which pre- 
vents the fringe pattern to vary with its maximum ampli- 
tude of 0.25 even in the case of optimum overlap, is a gen- 
uine effect of the dispersion-induced distortion between 
the early and the late wave packets. If the early and late 
wave packets were identical up to displacement, the visi- 
bility should not be affected, independently of their shape 
[23 |. It is mainly this criterion which decides whether a 
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single pulse momentum distribution is sufficiently well- 
behaved or not. Indeed, the joint detection probability 
(fTTjl can be evaluated analytically in the case of a Gaus- 
sian momentum distribution, revealing that the overall 
visibility is sensitive to the center-of-mass and relative 
motion momentum spreads (Tp,cm and crp,rci ^19]. Only 
if these are sufficiently small with respect to the relative 
momentum po, the resulting fringe pattern can exceed 
the visibility threshold 1/V2 as required for violating a 
Bell inequality. In our case, the momentum spreads are 
characterized by (Jp^tIpq = 0.024 and (Ap/po)'^ = 0.012, 
yielding a momentum distribution that is sufficiently well 
peaked to be capable of violating a Bell inequality, as can 
be seen in the lower plot in Fig. [3] We emphasize that 
the potential violation sustains over a sufficient number 
of fringes to be resolvable in experiment. 

Another viable approach to attain the expected fringe 
pattern as determined by a given (sufficiently well- 
behaved) momentum distribution is to fit it with a Gaus- 
sian. This permits to apply above mentioned analytic re- 
sult, as was done in [Toj . However, for spectra of the form 
(ITU)) such a least square fit has the tendency to slightly 
overestimate the accessible visibility. A more conserva- 
tive procedure than least square fitting would define the 
fit as an upper envelope of the corresponding momentum 
distribution. We note that the second momenta of the 
momentum distribution resulting from (jlOp are not use- 
ful for characterizing the widths of the peaks, while the 
Gaussian fitting procedure yields a comparatively good 
approximation in our context. 

A final remark should be made on the relative phase (jj-r 
between the early and late state component. As can be 
seen in pT|) . it also enters the joint detection probability 
and hence influences the fringe pattern. It was neglected 
in our numerical investigation since it only eff'ects a phase 
shift of the fringe pattern. Of course, this assumes that 
(j)T- remains constant from shot to shot, while an uncon- 
trolled variation of 0^ already on the order of lOOmrad 
would spoil the fringe pattern. In our Feshbach calcula- 
tion is determined by 

(j)r ^ [2UtT - fires ABT + Aires (i?rcs - Bo)t]/S + UgT. 

Hence, we have to require that the magnetic field pulses 
can be kept stable and are reproducible from shot to shot 



with a relative accuracy of 10 



V. CONCLUSIONS 

We described an experimental scenario for the gener- 
ation of a pair of dissociation-time entangled atoms and 
the subsequent interferometric Bell test. Assuming vi- 
able experimental conditions, the generated state of two 
counterpropagating atoms was shown to be capable of 
violating a Bell inequality. 

The experiment is based on the Feshbach-induced dis- 
sociation of a diatomic molecule, where the entanglement 
of the dissociating atoms is achieved by applying a se- 
quence of two magnetic field pulses. Hence, no compli- 
cated manipulation of the molecular state is required, 
the only experimental demand is to keep the molecules 
trapped in a controlled way and to carry out the magnetic 
field pulses with sufficient accuracy. The atom guides and 
the required switched, asymmetric Mach-Zehnder inter- 
ferometers are implemented by laser beams. Since only 
single-particle interference has to be performed, the in- 
terferometers on the opposite sides do not have to be 
kept stable with respect to each other. This allows to 
reach macroscopic distances between the two measure- 
ment devices, in our scenario on the order of 10 centime- 
ters. The procedure is completed by position measure- 
ments that only have to be able to distinguish the chosen 
beam splitter output port for each atom. This does not 
require prominent spatial and temporal resolution and is 
conveniently accomplished by resonant laser beams. 

The violation of a Bell inequality on the basis of 
dissociation-time entanglement would be the manifesta- 
tion of non-classicality in the motional state of two ma- 
terial particles, without resorting to internal degrees of 
freedom like spin. Here not only the two atoms are sepa- 
rated on a scale of several centimeters, but also the early 
and the late wave packets are separated by a distance of 
5mm on each side. Such a superposition of two consec- 
utively propagating wave packets for each atom, which 
would be established in this experiment, constitutes a 
macroscopic violation of the classical particle concept of 
localization and therefore presents a worthwhile feature 
of its own. 
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